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Driven Spin Systems as Quantum Thermodynamic Machines: Fundamental Limits 
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We show that coupled two level systems like qubits studied in quantum information can be used 
as a thermodynamic machine. At least three qubits or spins are necessary and arranged in a chain. 
The system is interfaced between two split baths and the working spin in the middle is externally 
driven. The machine performs Carnot-type cycles and is able to work as heat pump or engine 
depending on the temperature difference of the baths AT and the energy differences in the spin 
system AE. It can be shown that the efficiency is a function of AT and AE. 
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I. INTRODUCTION 

One of the main goals of thermodynamics has been to 
study heat engines and thermodynamic processes, dat- 
ing back to the now famous work of Sardi Carnot in 1824 
With the advent of nano physics and the control 
of quantum systems down to single atoms a better un- 
derstanding of thermodynamics on the basis of quantum 
mechanics is necessary. 

Since the first attempts to analyze thermodynamic ma- 
chines on the quantum level 0, Q considerable progress 
has been made in the last decades. Different kinds of 
models like machines built of harmonic oscillators, of un- 
coupled spins, particles in a potential or different three 
level systems have been studied [1, H, H, 0, Hi • 

Also the question about the violation of the second law 
of thermodynamics in the quantum regime has come up 
now and then. For a heat engine this would lead to an 
efficiency larger than the Carnot efficiency. All attempts 
to do so have failed and could be resolved, e.g., with the 
help of Maxwell's demon Q . 

Two level systems (TLS) like spins or qubits are the es- 
sential ingredients for quantum computation [10]. Much 
effort has been directed towards control of small clusters 
and chains of qubits in quantum optical systems 
nuclear magnetic resonance [12| and solid state systems 
[l3| . A serious problem in any such realization is the 
interaction of the respective quantum network with its 
environment. 

In the present work we study a model consisting of 
three TLS arranged in a chain in contact with two baths 
of different temperatures as studied for transport scenar- 
ios, e.g., in [ij. Ill, [l^. Here an energy gradient on the 
system and an incoherent driving of the TLS in the mid- 
dle let this system act as a thermodynamic machine. For 
possible experiments the set up may require more TLS's. 



Under special conditions the Carnot efficiency may be 
reached by a TLS heat engine but can never go beyond: 
If the Carnot efficiency is reached the machine flips its 
function, e.g., from a heat pump to a heat engine. 

We start with a discussion of the concept of work and 
heat. This is done by considering the change of the en- 
ergy expection value of a quantum system. With the 
help of the Gibbs relation heat can be associated with 
the change of occupation numbers of a quantum system 
whereas work is the change of the spectrum. 

We then introduce our thermodynamic machine con- 
sisting of three TLS's Thermodynamic properties 
can be imparted on this system by an appropriate em- 
bedding into a larger quantum environment [isl . fiol . [20l | , 
without the need of any thermal bath. In the present 
context, though, it is much simpler to settle for the open 
system approach based on a quantum master equation 
(QME). In Sec. lmBl the used QME will be introduced. 

Our numerical results are detailed in Sec. lIVI In Sec.lVl 
we compare the numerical investigation with an ideal 
TLS machine where ideal process steps are assumed. The 
obtained result is rather general and valid for any kind 
of TLS machines. 



II. THERMODYNAMIC VARIABLES 

A. Work and Heat 

To describe thermodynamic processes and machines 
one first has to define the pertinent variables heat, work, 
temperature and entropy for the system under consider- 
ation. Starting from the energy expectation value 



U ={E)= Tr {Hp} 



(1) 



for a quantum system H with discrete spectrum (p is the 
density operator) and considering the temporal change 

of {E) 
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(2) 
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change of work W can be associated with the first term 
of ^ where only the spectrum changes 



dt 



(3) 



W is the change per time of the i-th eigenvalue and is 
the corresponding occupation probability. The change of 
heat Q is then the second part of ^ 



dt 



Q-Tr 



(4) 



Equation ([2]) thus boils down to the famous Gibbs rela- 
tion 



AU = AW + AQ, 



(5) 






FIG. 1: Schematic representation of the system under investi- 
gation: An inhomogeneous 3-spin-chain is interfaced between 
two baths. Spin 1 (with energy splitting Si) and spin 3 (Ja) 
act as filters whereas spin 2 [(52(f)] is the working gas by de- 
formation of its spectrum. 



where AU is the energy change of the system. For cyclic 
processes work AW and heat AQ can also be calculated 
with the help of the S'T-diagram. For a closed path in 
the S'T-plane AU = and thus 

AW = -AQ = - jrdS. (6) 

While connected with bath a AQa can alternatively be 
calculated from the respective heat current Ja over one 
cycle of duration r 

AQ„ = / Jc^dt. (7) 
Jo 

Typically there are two baths a = h,c and thus two 
contributions (see Fig. [J) 

AQ = AQh + AQ,. (8) 



C. EfHciencies 

The efficiency of a heat pump is defined by the ratio 
of the heat AQh pumped per cycle to the hot reservoir 
and the work applied 

r,P = -AQh/ AW (11) 

which reduces for the Carnot heat pump to 

^Car - Th/in - Tc). (12) 

For the heat engine the efficiency is 

Tj^ = -AW/ AQh (13) 
which in the Carnot case leads to 

VL. = 1 - Tc/n. (14) 

III. DRIVEN SPIN SYSTEM 



B. Temperature and Entropy 



Hamilton-Model 



The temperature of a system can be defined if the state 
in the energy eigenbasis is canonical. For a TLS /x the 
temperature is given by 
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(9) 



with occupation probability of the energy level E'^. 
Due to the fact that all coherences will be damped out by 
the bath it is always possible to get a local temperature 
for a single TLS. The von-Neumann-entropy 

S^^-Tv {p,, In P,.} = -J2p'^ Inp; (10) 

i 

can then be taken as the thermodynamic entropy. 



The model under investigation is an inhomogeneous 
spin chain with nearest neighbor coupling of Heisenberg 
type described by the Hamiltonian 



E 



(15) 



The (T^'s are the Pauli-operators of the fith spin. A is the 
coupling strength which is chosen to be small compared 
to the local Zeeman splitting 6^, X 6^. Because 6^ ^ 
(5^+1 we call the spin chain inhomogeneous. 

We will need at least three spins in order to have this 
system work as a thermodynamic pump or machine. The 
spin chain is brought in contact locally with two baths at 
different temperatures as depicted in Fig. [TJ The detun- 
ing between spin 1 and spin 3 is (S13 = {Si — 5'i)/2 > 0. 
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B. Quantum Master Equation 

There are different ways to describe ttie tliermal behav- 
ior of quantum systems coupled to environments. Ex- 
amples are the path integral method [2l| or schemes 
based on the complete Schrodinger dynamics of the small 
system embedded into a larger quantum environment 
[isl . [igl . [20I I . Because it is much simpler for the present 
context we settle for a master equation approach. Such 
an approach has been widely applied to describe system 
bath models [22, 2||. 

To derive the master equation for our model one usu- 
ally starts from the Liouville-von-Neumann equation for 
the total system (we set h and the Boltzmann constant 
kg equal to 1) 



The Hamiltonian is composed of three terms 



H = H, + 



int 



(16) 



(17) 



with Hs the Hamiltonian of the relevant system, i/cnv 
the environment/bath Hamiltonian, Hi-at the interaction 
of coupling strength k between system and bath. With 
the help of a projection operator technique up to second 
order in k and with the use of the Born-Markov approxi- 
mation the dynamics of the reduced system density ps (t) 
reads 

d w ^ 
d^^^(^) ; 

- y dsTrenv |[i?int(i), [-^int (^ - s), p,{t) (g) Pcnv]]} , 

(18) 

where penv is a fixed state of the environment and Trenv 
denotes the trace over all degrees of freedom of the envi- 
ronment (see d^l)- 

In general the interaction Hamiltonian Hint is defined 

as 



Hi 



int 



>B,, 



(19) 



where Xi operates on the system and Bi on the envi- 
ronment. For the coupling with a single spin we take 
X = a^. By putting (fTO|) into (fT^ and going to the 
Schrodinger picture the following compact form can be 
obtained 



suppressing the system label s in the following. {l\ and 
m) are system eigenstates with the respective energy 
eigenvalue Ei/„^. ^{Ei — E„i) — ^{ujim) is the bath cor- 
relation tensor 



$(Wim) 



{B{s)B{0))ds 



(23) 



containing the bath correlation function 

(B(s)B(O)) - Trenv {B(s)B(0)pc„v}. (24) 
Assuming that the state of the bath is a thermal one 



Pc 



(25) 



(Zcnv being the partition function) and that the bath 
consists of uncoupled harmonic oscillators ^{ujim) takes 
the form 



^{uJhn) = K ( -^^^^ + e{uj,nl 



1 



(26) 



0{'^im) is the step function and the respective inverse 
bath temperature. 

For a three spin chain between two heat baths of dif- 
ferent temperatures Th and Tc and local coupling at the 
two chain boundaries with 



Xh = 0-1 ® I2 I3 
= il «) 12 ® 0-3. 

we get instead of ([20]) (cf. 0) 



-p = -\[H,p]+Vh{p)+Vrip). 



(27) 
(28) 



(29) 



The stationary state of (|20p for fixed 5^ is easily shown 

to be canonical of the form jf^'^'^ = e~'^^=/ Tr {e~'^^=}. 
However, with both baths in place the spin chain might 
be viewed as a molecular bridge generating a stationary 
leakage current — Jh ^ —Jc- Here the heat current J a 
between the 3-spin-system and the bath a can be defined 
by the energy dissipated via bath a (cf. [23 |) 



J„ = Tr{i?P„(p)}. 



(30) 



In the following a current out of bath a into the machine 
will be defined as positive. 



(20) 



As in [1^ we use the dissipator Vips) 

Vip)^[X,Rp] + [X,Rp]^ (21) 

with 

(l\R\m} ^ {l\X\m)^Ei-E^) (22) 



IV. NUMERICAL RESULTS 



A. Non-equilibrium stationary states of a spin 
chain 



First we note that the heat current through a spin 
chain depends on the local Zeeman splittings within the 
system. To analyze the heat current we solve (|29l) and 
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FIG. 2: Stationary heat current Jh [see (I30|l ] (from hot bath) 
and Jc (to cold bath) as function of the local energy splitting 
52 of spin 2 with 5i = &3 = 1, Th = 2.63 and = 2.5. 
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FIG. 3: Heat currents Ja{t) for the heat pump over one cycle 
with duration t = 2-k/lo = 804.25 for Tu = 2.63 and Tc = 2.5. 
The peaks are a result of the resonance-effect. 



calculate the stationary state of the system p'^'^* . With 
this solution and with the help of (|30p we can then cal- 
culate the currents Jq for each bath. 

We consider a system with 5i = 5^ ~ 1. Both heat 
currents (|5D|) as function of 62 are shown in Fig. O Jh is 
positive and the relation Jh — ~Jc is fuUfilled. If 82 — 
81 = (^3 (the homogeneous case) the heat currents reach 
their maximum. By detuning the local energy splitting it 
is thus possible to uncouple the respective bath from the 
rest of the system. This resonance effect will now be used 
to build out of three spins a quantum thermodynamic 
machine. 



B. Time dependent behavior: spin system as 
thermodynamic heat pump or machine 

1. The heat current 

The above situation changes when the energy splitting 
of spin 2 is chosen to be time-dependent, i.e., 

52{t)^sin{ujt)h + hfi. (31) 

Different energy splittings of the boundary spins, e.g., 
5i = 2.25 and 63 = 1.75, are used to install left/right- 
selective resonance effects. The parameters of (|3T|) are 
chosen as bo = {Si + 5z)/'2. and h is the detuning ^13. To 
enable the bath to damp the system lj <^ 82 must be 
fuUfiUcd. 

For solving we have used a four step Runge-Kutta- 
algorithm. At each time step the bath correlation func- 
tion is calculated explicitly. We choose the following pa- 
rameters for our numerical results: A = 0.01, k = 0.001, 
5i = 2.25, (53 = 1.75, uj = 1/128, = 2.5 and Th is var- 
ied, unless stated otherwise. Both coupling parameters 
A and n are chosen to stay in the weak coupling limit. 

Now when spin 2 is driven periodically as in pip , we 
can distinguish four different steps: 



1. Spin 2 (the "working gas") is in resonance with 
spin 3 [52 (t) ~ 5^] and thus couples with bath c at 
temperature T^- Because of this energy resonance 
the current Jc via spin 3 will be large, whereas the 
current Jh via spin 2 will be negligible. The occu- 
pation probabilities of spin 2 and 3 approach each 
other and so do the respective local temperatures. 

2. Quasi-adiabatic step: Spin 2 is out of resonance 
with spin 3 [5i > S2{t) > S3], now J^, is suppressed 
while Jh nearly stays unchanged. The occupation 
probability of spin 2 does not change significantly 
and there is almost no change in the entropy 82- 

3. Spin 2 is in resonance with spin 1 [62(1) « Si] and 
by that in contact with bath h at temperature Th- 
Jh is large whereas Jc is very small. The local tem- 
peratures of spin 1 and 2 nearly equal each other. 

4. Quasi-adiabatic step, as in step 2. 

Fig. [3] shows the heat currents Jq of both baths over 
one period with the bath temperatures Th = 2.63 and 
Tc = 2.5. As can be seen the resonance effect decouples 
spin 2 from the bath if its energy splitting is different 
from the boundary or filter spins. This decoupling is 
never perfect, though. As a consequence there is a leak- 
age current Jl which will be discussed in more detail 
later on. 



2. Heat, Work and Efficiencies 

That the studied system indeed works as a heat pump 
can be seen from the S'2T'2-diagram of spin 2 in Fig. U) 
The local entropy S2 of spin 2 is given by pH)) and the 
local temperature T2 by @. The four different steps as 
explained in Sec. lIVBTl are shown as well as the direction 
of circulation. 

To determine the efficiency of this heat pump one needs 
to know the quantity of heat AQh pumped to the hot 



T2 

FIG. 4: S'2r2-Diagram for the quantum heat pump for Th = 
2.63, Tc = 2.5 (AT = 0.13) and r = 27r/tj = 804.25. The 
arrows indicate the direction of circulation. 
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FIG. 5: S'2r2-diagram for the quasi-statically driven quantum 
heat pump (with parameters as in Fig. |4]). Because of the 
leakage current the enclosed S'2r2-area vanishes and no work 
is exchanged. 

bath and the used work ISW . AQh can be calculated by 
integrating the heat current Jh over one period [cf. ([7])]. 

The exchanged work AVF is given by the area enclosed 
in the 5'2T2-plane according to ([6]). We find that indeed 
AW+AQc+AQh = in all cases and by that confirm the 
use of T2 and 5*2 as effective thermodynamic variables. 

In contrast to the Carnot model our machine is work- 
ing in finite time. If driven too fast the bath is not able to 
damp the system and if driven too slow (quasi-stationary) 
the system would have reached its momentary steady 
state transport configuration (i.e. a; ^ k). The S2T2- 
area then vanishes as depicted in Fig. [51 This is caused 
by the leakage current. 

Figure [7] shows the Carnot efficiency for the heat pump 
^Car' ^^'^ machine rj^^^ and the respective efficiencies 
for our quantum heat pump ry^j^^ [according to pTjl ] and 
machine rj'^^ [according to as a function of the tem- 
perature difference AT = Th — T^- We point out the 
following interesting findings: 

• The efficiency curve of the quantum heat pump or 
machine is always below the respective Carnot effi- 



T2 

FIG. 6: S'2T'2-Diagram for the quantum heat engine (AT = 
0.83 > ATmax with Th = 3.33, = 2.5 and r = 2-K/ijj = 
804.25. The arrows indicate the direction of circulation. 
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FIG. 7: The Carnot-efBciency rj^^^ and the efficiency r]^^ 
of the quantum heat pump (AT < ATmax) and tjq^^ and 
77qjn of the heat engine ( AT > ATmax) as function of the 
temperature difference AT. Following parameters are chosen: 
Tc = 2.5, 5i = 2.25, ^3 = 1.75 and r = 27r/tj = 804.25. 



ciency. As expected, the 2nd law is never violated. 

• For AT = 77^,-,-^ does neither diverge nor go to 
zero. This means that the machine can start out of 
equilibrium and begin to cool a reservoir. 

• At a specific temperature difference AT, here 
ATniax ~ 0.6, the heat pump switches to operate as 
a heat engine. To illustrate this fact Fig. [S] shows 
the area in the S'2T2-plane for AT = 3.33 > AT^ax- 
As depicted, the direction of circulation has re- 
versed. 

To make the last point more plausible Fig. [8] shows the 
work AW , the heat Qh and Qc as function of AT. While 
AT is increasing, AQh and AQc are decreasing as well 
as AW until first AQc changes its sign, then Qh and last 
AW. At the point where AW = (for AT = AT^ax) 
only the leakage current Jl is flowing from the hot bath 
to the cold one. Beyond this ATmax the system starts to 
work as an engine. 
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FIG. 8: Heat AQc and AQ^ and work AW performed over 
one cycle as function of the temperature difference AT (same 
parameters as in Fig. [T] The inset shows these functions 
around the point AT = ATmax in more detail. AQl is the 
leakage heat per cycle. 



In contact with spin 1 spin 2 exchanges heat AQjf with 
the hot bath at T^. No work will be done and only the 
occupation probabilities of spin 2 will change to a ther- 
mal state with T2 = T^. The exchanged heat can be 
calculated by the energy difference before and after ther- 
malisation 



5i 



tanh 



2Tc 



tanh 



2T,, 



(34) 



Then spin 2 is driven back to the energy splitting of spin 3 
[E^ ^ El = ^ and El ^ El = The work VFi^g for 
this step is given by 



^.(.,-.3) tanh (A.) 



Finally the heat Q^^ 

-53 



AQ^^ = 



— tanh 



2Tr 



(35) 



(36) 



V. ANALYTICAL RESULTS 

A. Ideal Quantum Machine 

To understand the above numerical results wc com- 
pare them to the maximal reachable heat and work which 
could be pumped or extracted by a TLS quantum ma- 
chine. All process steps will be taken to be ideal steps. 
By ideal we mean that we have total control of each pro- 
cess step. Then no leakage current will disturb the system 
and the heat exchange at bath contact will be without 
loss. 

In addition we assume a machine which only works 
during the adiabatic steps. Heat will be exchanged only 
if the machine is in contact with a bath. This can be 
compared with the Otto-cycle d, HI] . 

We start with spin 2 in contact with spin 3 and thus 
with cold bath. The state of spin 2 after this contact is 
a canonical one of the form 



1 /e*3/(2T.) 
Z \ e-^3/(2T.) 



(32) 



Z is the partition function and we have assumed that the 



El 



and the 



energy of the ground state is E2 
excited state E\ = E^ = ^ because both spins are in 
resonance. 

After this equilibration with the cold bath at Tc the 
spin 2 is driven until its local energy splitting is equal to 
spin 1 (i;^ = i;? ^ and El ^ El = ^-). The work 
for this step can be calculated with ([3]). This step is 
adiabatic as p2 does not change. The work W3_>i is then 
given by the energy difference before and after reaching 
the splitting of spin 1 



l^3^i = ^('53-'5i)tanh('A 



(33) 



will be exchanged with the cold bath via spin 3. 
total work AWtot is given by 



Al^tot = W^3-.l 

The Gibbs-relation 



The 



(37) 



(38) 



can easily be verified. 

With the help (|33l) - ([37]) it is now possible to calcu- 
lated the efficiency of this ideal machine. For the heat 
pump we get 



id 



V 



for the machine 



AQk' 



5i - 83 



5i - S3 



Si 



(39) 



(40) 



This result is similar to that obtained by Kieu [2^ and 
is the maximum a TLS can reach. Here we want to com- 
pare the efficiency of the ideal pump rj?^ and engine 77?) 
with the respective Carnot efficiencies for the parameters 
used for our numerical results. 

Figure [5] shows the Carnot efficiencies as well as the 
one from ([55]) and ([3D]) . rifj^"^ is always below r]^^^^ until it 
reaches a maximal temperature difference ATmax (with 

= 2.5, Si = 2.25 and (53 = 1.75 we get AT„iax = 0.714). 
At this temperature the heat pump is working loss- 
less and no heat can be pumped. Just like the quasi- 
stationary Carnot heat pump this pump has zero power. 
Only in this particular case rj^^ = rjQ^j.. By further in- 
creasing the temperature Th the heat pump starts work- 
ing as a heat engine. 

Figure [TUl illustrates this behavior where AM^"*, AQ'^ 
and AQjf are depicted as function of AT. At ATmax no 
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FIG. 9: Carnot-efficiency rj^^^ for the heat pump and engine 
'7car as function of temperature difference AT while Tc — 
2.5, 5i = 2.25 and ^3 — 1.75 as in Fig. [T] -q^^-i and rj^^ are 
the efficiencies of the ideal pump/engine [see (|39} and (|40|) ]. 
?7?j = 12.36 and AT^ax = 0.22 can be reahzed for 61 = 1.904, 
53 = 1.75 and Tc = 2.5. 



FIG. 11: Fitted efficiency rj^^ for the quantum heat pump and 
quantum engine rj^^ as function of temperature difference AT 
while Tc = 2.5, 5i = 2.25 and Ss = 1.75 as in Fig. [7] vf^ and 
rj^^ are the efficiencies of the ideal pump / engine [see (|39} and 
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FIG. 10: Work AW'"^, heat AQ',^ from/to the hot bath and 
heat AQjf from/to the cold bath for the ideal machine as 
function of temperature difference AT while Tc — 2.5, 5i — 
2.25 and ^3 = 1.75 as in Fig. El At AT = ATn^ax < = Veil 
and therefore AW''^ = 0, AQ'^'^ = and AQ'^ = 0. 



heat ^Q]f is pumped and therefore no work used or no 
heat exhausted to do work. 

This is quahtatively the same behavior as our model 
shows in Fig. [7] and Fig. [51 Two differences can be seen. 
First the critical temperature in our numerical result de- 
viates from the theoretical expected one. From the nu- 
merics we get ATmax ~ 0.6. Second the inset in Fig. [5] 
shows that AQc changes its sign before AQh does. The 
reason for both effects is due to the leakage current as 
will be explained below. 

For a given bath temperature (like in our example Tc = 
2.5) it is possible by changing the energy splittings of 61 
and/or S3 to influence ATmax- In Fig. [9] also a different 
efficiency fjf^ is depicted. 77?^ can be realized by increasing 
i5i so that ATmax will be decreased to ATmax- 



B. Quantum machine with leakage current 

The efficiency of an ideal two level quantum machine 
is independent of AT except at AT = ATmax, where it 
jumps between its heat pump and its heat engine value. 
The efficiency obtained from the numerical simulation 
deviates somewhat from this expected behavior. For the 
heat pump the efficiency of our model is even larger than 
the ideal one (see Fig. [TT|) . To understand this effect 
we analyze the leakage current from a phenomenological 
point of view. 

First we assume that the leakage current causes the 
gas spin 2 to approach a thermal state which is not in 
accordance with the bath temperature. In this case AQh 
and AQc will be decreased. This effect is responsible for 
the vanishing of ?7q4™ before reaching ATmax- But it can 
not explain why the efficiency of rjP^ is sometimes larger 
than 77P . 

Taking into account that also less work is performed 
due to the leakage current it is possible to find a larger 
efficiency. This can be interpreted in that the gas spin 2 
does not "see" the full energy splitting 61. As shown in 
Fig. [TT] our phenomenological model fits the numerical 
data quite good. 

For the efficiency of the heat engine rj^^^ it can be seen 
from Fig. [TT] that it is always worse than the ideal engine 

V ^ e 



VI. CONCLUSION 

We have studied a driven 3-spin system coupled to two 
split heat baths. We have shown that such small quan- 
tum networks may be used not only as quantum infor- 
mation processors but also as quantum thermodynamic 
machines. For the latter proposal we would primarily ex- 
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ploit the (time-dependent) deformation of discrete spec- 
tra and associated resonance transfer. 

While interesting functionahty appears already for N 
= 3 spins, also larger spin networks subject to such very 
limited control could be envisaged without losing inher- 
ent stability: eventually this stability is dictated by the 
increase of entropy, i.e., by the second law of thermody- 
namics. 

For a thermodynamic TLS machine working with ideal 
heat transport and adiabatic steps we have derived an 
ideal efficiency. This efficiency is independent of the bath 
temperatures. By tuning the energy splitting of the TLS 
the quantum thermodynamic machine can be used as a 
heat pump or heat engine. The Carnot efficiency will only 
be reached when a TLS machine is working losslessly. 

Taking dissipation into account it is possible to under- 
stand the leakage current present in our numerics from a 
phenomenological point of view. Surprisingly a leakage 
current could even increase the efficiency of a heat pump 
whereas for a heat engine it only decreases the efficiency. 

There are a number of different options for implemen- 
tations [1^ [l^l and also various possibilities to introduce 
the time-dependent control. For simplicity we have re- 
stricted ourselves here to external driving; alternatively 
one might look for autonomous system design s I28l| . e.g., 
by using a mechanical oscillator (cantilever) [23]. Arti- 
ficial autonomous nanomotors powered by visible light 
have recently been demonstrated experimentally [3(1 |. 

From a fundamental point of view several interesting 
questions remain: What is the status of thermodynamic 



variables for such quantum systems? To what extent are 
these measurable in the nano-domain - without being op- 
erators? And if measured, how would the measurement 
result fluctuate j3f|? 

As noted already, a two-level system diagonal in its 
local energy basis can always be described as canonical 
with some temperature T, i.e., there is conceptionally no 
space for nonequilibrium here. It is remarkable that for 
periodic operation work can then be associated with the 
area defined by the closed path in the effective entropy- 
temperature plane for the driven spin, as in macroscopic 
models. 

This may challenge the subjective ignorance interpre- 
tation of none-pure states as classical mixtures, i.e., as- 
suming the individual spin to be either up or down at any 
time. If the thermal state was taken to result from quan- 
tum entanglement with the environment ,18;], this clas- 
sical picture would no longer be needed; those concepts 
from quantum information seem to be more appropriate 
here. 
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